Subdiffusion on a Fractal Comb 
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Subdiffusion on a fractal comb is considered. A mechanism of subdiffusion with a transport 
exponent different from 1/2 is suggested. It is shown that the transport exponent is determined by 
the fractal geometry of the comb. 
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A comb model was introduced for understanding of 
anomalous transport in percolating clusters 0, 0| and 
it was considered as a toy model for a porous medium 
used for exploration of low dimensional percolation clus- 
ters 0, 3, as well. It is a particular example of a 
non-Markovian phenomenon, which was explained in the 
framework of continuous time random walks 
the last decade the comb model has been extensively 
studied for understanding of different realizations of non- 
Markovian random walks both continuous and dis- 
crete fill ]. 

Usually, anomalous diffusion on the comb is described 
by the 2D distribution function P = P(x, y, t), and a spe- 
cial behavior is that the displacement in the ^-direction 
is possible only along the structure axis (x-axis at y = 0). 
Therefore, diffusion in the x-direction is highly inhomoge- 
neous. Namely, the diffusion coefficient is D xx = DS(y), 
while the diffusion coefficient in the y-direction (along 
fingers) is a constant D yy = D. Therefore, this inho- 
mogeneous diffusion is described by the Fokker-Planck 
equation in the dimensionless time and coordinates 



L FP P(x,y,t) = d t P - 5{y)dlP 



d 2 P = 0. 



(1) 



It is obtained by the rescaling with relevant combina- 
tions of the comb parameters D and D, such that the 
dimensionless time and coordinates are D 3 t/D 2 — > t 
Dx/D — ¥ x, Dy/D — > y, correspondingly [12j . 

The fractional transport along the structure x axis is 
described by the transporting contaminant distribution 
p(x, t) — P(x, y, t)dy. It was shown that Eq. ([T]) 
is equivalent to the fractional Fokker-Planck equation 



dfp(x,t)- ~d 2 p(x,t) = 0, 



(2) 



from where subdiffusion can be immediately obtained: 

J x 2 p(x, t)dx ~ y/t. Here 9 t 2 is a fractional time deriva- 
tive, which is a formal notation of an integral with a 
power law memory kernel. For < a < 1 it reads 



d?f(t) 



T(l - a) 



-d r f(r)dr. 



(3) 



Subdiffusive mechanism with an arbitrary transport ex- 
ponent was also suggested by changing either the bound- 
ary conditions for diffusion in the fingers 13 16], or in- 
troducing a dependence of the diffusion coefficient on 



time and space [l7| . In this paper we consider a frac- 
tal comb, when diffusion is highly inhomogeneous along 
the y fingers, as well. Namely, it takes palace for those 
coordinates of the x axis, which belong to a fractal set 
S„(x), and is defined by a characteristic function x(x), 
such that Dyy = D\(x), where x( x ) = 1, if £ € S v (x) 
and x( x ) — 0, if x ^ S u {x). The fractal set S v {x) is a 
random fractal with a fractal dimension < v < 1 em- 
bedded in the ID Euclidian space (of the x axis). Such 
generalization of the comb model to a discrete (fractal) 
comb model for consideration of fractional transport in 
discrete systems is more realistic situation for theoretical 
studies of transport properties in discrete systems with 
complicated topology including fractal ones like porous 
discrete media [l8j . electronic transport in semiconduc- 
tors with a discrete distribution of traps, cancer devel- 
opment with definitely fractal structure of the spreading 
front , see e.g., reviews @,[Hj], and infiltration of diffusing 
particles from one material to another (20j . 

Hence, we study the following dimensionless equation 



dtV - 5{y)d 2 x V - X (x)dyV = . 



(4) 



The initial condition is V(x,y,0) = 6(x)5(y), and 
the boundary conditions on infinities have the form 
■p(±oo, ±oo, t) = and the same for the first deriva- 
tives with respect to x and y V x (±oo,±oo,t) — 
V y (±oo,±oo,t)=0. 

Our main purpose is to evaluate the second moment 



(x 2 )(t)= / x 2 V(x,y,t)dxdy 



(5) 



as a function of time. Therefore, the forthcoming anal- 
ysis of Eq. (g]) is supposed to be carried out under the 
integration sign. Using properties of the characteristic 
function 

X 2 ( X )=X( X ) and d x x(x) = , (6) 

a solution of Eq. (0| can be presented in the form 

V{x,y,t) = X {x)P{x,y,t), (7) 

where P(x, y, t) is a solution of the continuous comb 
model and we shall show that an equation for this func- 
tion coincides with Eq. ([1]). But first, one should under- 
stand a physical meaning of the distribution V(x,y,t) 
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FIG. 1: Comb structures. The upper strip corresponds to 
the continuous comb model. The lower strip is a sketch of 
the fractal comb with a specific distribution of fingers, which 
corresponds to the one third Cantor set (at the forth step 
of the construction). Asking a question: "one third of what 
scale on the infinite line is it?", one should recognize that 
this construction is not representative. Therefore, the fractal 
set F u can be considered as a random fractal distribution 
of the fingers without any specifically defined construction 
algorithm. It is worth admitting that the fractal distribution 
of teeth, in general case, is multifractal, and the uncertainty 
of this construction should be stressed as well. 



based on properties of the characteristic function x( x )- 
While the first property in Eq. ([5]) is obvious and follows 
from the definition of x( x ): the second expression de- 
serves an explanation. To show this, let us consider the 
iVth step of the fractal set S v construction. It is a union 
of disjoint intervals Axn- In general case of a random 
fractal, these are random intervals. In the limiting case 
one obtains S v = limAr_ i . 00 (J Ax^. Therefore, the char- 
acteristic function on every interval Ax j — [xj , x j + Aa: at] 



Ax 



N 



Differentia- 



is x(AxAr) = 8(x — Xj) — Q(x — 

tion of the characteristic function on every interval yields 
J^x(A- x n) — S(x — Xj) — 5(x — Xj — Axn)- In the limit 
N — > oo it tends to zero (under the integration), since 
P(x, y, t) and its derivatives are continuous functions. 

Now we can return to the continuity property of 
V(x, y, t) that can be understood from the calculation of 
the second moment (x 2 )(t) in Eq. The presence of 
the characteristic function in this expression means that 
the integration is performed over the fractal volume [2l| . 



It means that J x{ x )dx —> f \x\ u 1 dx ~ Wy+l) \ x \"> 
where T(y) is the Gamma function. This yields the ex- 
pression for (x 2 )(t) 



(x 2 )(t) 



1 



I» J- 



x 2 \x\"- 1 P(x,y,t)dxdy . (8) 



Using this smoothing procedure, we can show that an 
equation for the distribution function P(x,y,t) is Eq. 
(HJ). We have from Eqs. (g]) and properties ^ of the 
characteristic function x( x ) that for any arbitrary func- 
tion f(x) 



\x\»- 1 f{x)L FP P(x,y,t)dx = Q. 



(9) 



Therefore, we have equation LppP(x,y,t) = 0, which 
exactly coincides with Eq. ([T]) and is valid for all x. 

Now we at a position to determine (x 2 )(t). Taking into 
account Eq. ©, one obtains from integration of Eq. © 
over y with f(x) = x 2 



d t (x 2 )(t) = 



I» 



\x\ u - 1 x 2 dlP(x,0,t)dx. (10) 



Here we take into account that d 2 P(x,y,t) = 
due to the boundary conditions. A relation between 
P(x, 0, t) and P(x, y, t) can be established in the Laplace 
domain. Performing the Laplace transform CP{x, y, t) = 
P(x, y, s) in Eq. (jTJ) , it is readily to see that P(x, y, s) — 
P(x, 0, s)e~^\ y \ satisfies the equation. After integrating 
over y, it yields 



P(x,0,s) 



P(x,V, s)dy= -Vsp(x, s).. (11) 



This result can be taken into account after the Laplace 
transform in Eq. (|10[) . that yields an expression for 
the second moment in the Laplace domain C[(x 2 )(t)} — 
(x 2 )(s). It reads 



s(x 2 )(s) 



I» J- 



x\ 1+v d 2 p{x, s)dx , 



(12) 



where p(x, s) = ^^ 3/2 exp(— V2s 1 / 2 |a;|) can be obtained 
from Eq. ([2]). After the Laplace inversion one obtains 
the second moment 



(x 2 )(t) =K v t^ 



(13) 



")/r(* 



, , | )T{v)y/2 1+v is a generalized 
Finally, we obtain subdiffusion on 
t M with the transport exponent \ < 
= 1, one observes subdiffusion with 



where K v = T(2 
diffusion coefficient 
the comb (x 2 )(t) 
(i < \. When ; 
fj, = 1/2. To obtain subdiffusion with i < [i < 1, one 
considers advection along the structure x axis instead of 
diffusion. This yields for the transport exponent to be 

A* 2 • 
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The main deficiency of the obtained result in Eq. (jT3|) 
is that it is based on the presentation of the probability 
distribution function as a product of a continuous func- 
tion and the characteristic function in Eq. (JT]) . Although 
the inferring of Eq. (fT3|) is correct, this presentation 
can leads to wrong result, because the probability dis- 
tribution function V{x, y, t) must be continuous at every 
point. To overcome this deficiency, we refuse the locality 
property. To this end, the following procedure of coarse 
graining of the Fokker-Planck equation ^ is suggested. 
First, we apply the Fourier transform to Eq. §Q with re- 
spect to the x coordinate. To apply this transformation 
to the last term in Eq. (j4]) , we use the following auxiliary 
identity 

x(x)f(x) = d x ( x (y)f(y)dy , x(o) = o . 



Here for brevity we define f(x) = V(x,y,t). This in- 
tegration with the characteristic function can be carried 
out by means of a convolution [23 |. Note that 



x(y)f(y)dy 



E 

Xj£S„ [0,x] 



f(y)S(y-x j )dy, (14) 



where we use that 



E % 



A*) 



x i £S v 



is a fractal density, such that on the finite interval (0, x), 
the integral J Q dfj,(y) ~ x v corresponds to the fractal 
volume. Therefore, Eq. (|T4l reads 

x(y)f{y)dy = / f{y)dn{y) . 
o Jo 

The last expression can be rewritten as a convolution 
integral. Due to Theorem 3.1 in Ref. [23[, we have 

f(y)dfi(y) ~ f\x - yT' 1 f{y)dy , (15) 

o r(i/) 7 

where A v is a constant, defined by the conditions of the 
theorem. In sequel we disregard this parameter, putting 
A v = 1. This integration is a Riemann-Liouvillc integral 
(see, e.g., gUl) 

Wf{x) = ^dx J\x-yT- 1 f{y)dy. 

Here we use a standard notation Ql^f(x) to define in- 
tegration with a power law kernel, see Eq. ([3]). The 
integration can be presented in the form of the inverse 
Laplace transform Cf(x) = f(s), which reads 

„£/(*) = £- l £[oiZf{x)\ = c-h-'m . 

Therefore, after the variable change s — iz, the Fourier 
transform of the last term in Eq. Q yields 



where T x f{x) = P(k,y,t) — f(ik). One takes into ac- 
count that the result should be symmetrical with respect 
to the negative x < 0. Therefore, the Fourier transform 
of Eq. flU) yields 

d t V(k, y, t) = -5{y)k 2 V{k, y, t) + ^ d 2 y V(k, y, t) . 

(17) 

Now we perform the Laplace transform with respect to 
time CV{k, y, t) = P(k, y, s) = G(k, y, s). This yields 

sG = -5{y)k 2 G+\k\ 1 - l/ d 2 y G + 8{y) (18) 

with the solution 



G(k,y,s) = eM-\y\Vs\k\»- 1 )g(k,s), (19) 

where g(k,s) = G(fc,0, s) is the Fourier-Laplace image 
on the structure axis at y = 0. Again, we are interest- 
ing in dynamics along the structure axis by studying the 
probability distribution function V\ (see Eq. (jTTJ) ) : 



Vi{x,t) 



V(x,y,t)dy. 



Integrating Eq. (|19|) over y one obtains 



G(fc,0,s) 



G(k,y,s)dy. 



(20) 



(21) 



Therefore, integrating Eq. (fT5|) over y yields the 
Montrall- Weiss equation that, after the Fourier and the 
Laplace inversions, reduces to the fractional Fokker- 
Planck equation. It is a particular case of a general equa- 
tion 

a?Pi{x,t) = ±V%P 1 (x,t), 0<a<l, (22) 

where /3 = | + S. It describes a competition between 
long rests and long flights. We stress that when a = 1/2, 
it corresponds to the comb model, see Eq. ©. Here we 
use the formal definition for the Riesz-Weyl fractional 
space derivative in form of the Fourier inversion (see, 
e.g., SE1): 



Vlf{x)=f- l [\k\»f{k)]. 



(23) 



This equation was studied in [23, |27( (see also review 
@). The correct form of the mean squared displace- 
ment, which estimates the competition of "laminar mo- 
tion events" (flights) and "localization" (waiting) events 
in the Levy walk picture was obtained through the rela- 
tion [HI valid for Eq. ([22]) 



(24) 



(x 2 (t)} ~f = i x 



F x [ X {x)V{x,y,t)} = {ikf- v V{k,y,t) 



(16) 



One easily checks that this result with fj, = 1 — u/2 has 
a correct limit for v = 1, when it corresponds to the 
continuous comb model with /i = 1/2. 

In conclusion, we presented two approaches to study 
subdiffusion on the fractal comb, when the fractal trap 
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distribution is determined by the characteristic function 
x(x). In this case, it is tempting to look for a solu- 
tion in the multiplicative form of Eq. 0. Thus the 
forthcoming analysis, based on this presentation of the 
probability distribution function is rigorous. The main 
deficiency of this approach is that it violates the continu- 
ity property of the probability distribution function. To 
overcome this deficiency, a coarse graining procedure of 
the Fokker-Planck equation ((4]) is suggested. It is based 
on the possibility to perform the Fourier transform for 
Eq. ((3]) exactly. We obtained that inhomogeneous, frac- 
tal distribution of traps in the comb model leads to Levy 
jumps that complicates fractal diffusion and leads to the 
competition between long jumps and localization inside 
traps. This phenomenon is described by the fractional 
Fokker-Planck equation (|2"2"j) . As a result of this compe- 
tition, subdiffusion, which is the dominant process, real- 
izes. 



We admit that the results of the either approach of 
Eqs. (TTBI or (f2"4"|) have correct limits for v = 1. A spe- 
cific property of Eq. (|24|) is that for v = 0, it corre- 
sponds to normal diffusion with fi = 1. While the first 
limit with v = 1 is well understood, the second one is 
not so obvious. Indeed, when the Hausdorff dimension is 
v = 0, there are no traps, and normal diffusion is antici- 
pated. Nevertheless, this behavior results from the frac- 
tional Fokker-Planck equation (|22|) with 3/2 fractional 
space derivative and 1/2 fractional time derivative. This 
is a special point of a transition from subdiffusion to su- 
perdiffusion 0, [26|. For the comb model subdiffusion 
is the dominant process, there is no superdiffusion, and 
jj, = 1 is the boundary point. 
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